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O ■ ABSTRACT 
CO 

We investigate radial and vertical metallicity gradients for a sample of red 
| clump stars from the RAdial Velocity Experiment (RAVE) Data Release 3. 

We select a total of 6781 stars, using a selection of colour, surface gravity 
and uncertainty in the derived space motion, and calculate for each star 
' a probabilistic (kinematic) population assignment to a thin or thick disc 

using space motion and additionally another (dynamical) assignment using 
stellar vertical orbital eccentricity. We derive almost equal metallicity gra- 
dients as a function of Galactocentric distance for the high probability thin 
disc stars and for stars with vertical orbital eccentricities consistent with be- 
ing dynamically young, e v ^ 0.07, i.e. d[M/H]/dR m = -0.041 ± 0.003 and 
dlM/Hj/dRm = -0.041±0.007 dex kpc" 1 . Metallicity gradients as a function 
of distance from the Galactic plane for the same populations are steeper, i.e. 
d[M/H]/dz max = -0.109 ± 0.008 and d[M/H]/dz max = -0.260 ± 0.031 dex 
kpc -1 , respectively. R m and z max are the arithmetic mean of the perigalactic 
and apogalactic distances, and the maximum distance to the Galactic plane, 
respectively. Samples including more thick disc red clump giant stars show 
systematically shallower abundance gradients. These findings can be used to 
distinguish between different formation scenarios of the thick and thin discs. 
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1 INTRODUCTION 



Metallicity gradients play an important role in understanding the formation of disc pop- 
ulations of the galaxies. In the Milky Way Galaxy, there is extensive information estab- 



lishing a radial gradient in young stars and in the interstella r medium ( jShaver et al. 



1983 



Luck. Kovtyukh fc Andrievsky 



2006 



Luck fc Lambert 



201 ll ) . Values typically derived are 



d[F e / H] / dRc = —0.06 ± 0.01 dex kpc -1 , within 2-3 kpc of the Sun. Much effort to search 
for local abundance variations, and ab undance variatio ns with azimuth, show little if any 
detectable variation in young systems ( ILuck et al.ll201l[ ). showing the inter-stellar medium 
to be well-mixed on quite large scales, supporting the importance of gas flows. 

Extant data suggest vertical metallicity gradients in the —0.4 < d[M/H]/dz < —0.2 



dex kpc 1 range for relatively small distances from the Ga 



(Trc 



: zger. Pel fc Gabi 



1995 



Karaali et al. 



2003 



Du et al. 



2004 



actic plane, i.e . z < 4 kpc 



Ak et al. 



2007a ; 



Peng. Du & Wu 



201 ll ). For intermediate z distances, where the thick disc is dominant, the vertical metallicity 
gradient is low, d[M/H]/dz = —0.0 7 dex kpc" 1 and the radial gradient is only marginal, 
-0.02 ^ d[M/H]/dR ^ dex kpc" 1 Jllong et al.lhoOlf ). There is some evidence that metal- 
licity gradients for relatively short vertical distances, z < 2.5 kpc, show systematic fluc- 
tuations with Galactic longitude, similar to th ose of thick-disc scaleheight, which may be 



interpreted as a common flare effect of the disc ( 



Ak et al 



2007b 



Bilir et al. 



2008 



Bilir et al 



Yaz&Karaal 



20101). 



2006 



Cabrera-Lavers et al. 



2007 



Quantifying the abundance distribution functions and their radial and vertical gradi- 
ents in both thin and thick discs can be achieved using stellar a bundances, especially t hose 
from major surveys such as RAdial Velocity Experiment ( RAVE; 



Sloan Digital Sky Survey (SDSS; 



Abazajian et al. 



Steinmetz et al. 



2006|) and 



20031 ). RAVE is a multi- fibre spectro- 



scopic astronomical survey in the Milky Way, which covers just over the half of the Southern 
hemisp here, using the 1.2-m UK Schmidt Telescope of the Australia n Astronomical Obser- 



vatory (iSteinmetz et al. 



2006 



Zwitter et al 



2008 



Siebert et al 



201l[ ). RAVE's primary aim 



is to derive the radial velocity of stars from the observed spectra for solar neighbourhood 



* E-mail: sbilir@istanbul.cdu.tr 
f Retired 



Metallicity Gradient from RAVE DR3 3 

stars. Additional information is also derived, such as photometric parallax and stellar at- 
mospheric parameters, i.e. effective temperature, surface gravity, metallicity and elemental 
abundance data. This information is important in calculating metallicity gradients, which 
provides data about the formation and evolution of the Galaxy. As the data were obtained 
from solar neighbourhood stars, we have limitations to distance and range of metallicity. 
However, the metallicity measurements are of high internal accuracy which is an advantage 
for our work. 

In a recent study carried out by 



Chen et al. 



(120111 ). a vertical metallicit y gradient of 



-0.22 dex kpc 1 was claimed for the thick disc. They used the SDSS DR8 (lAihara et al 



20 111 ) data set to identify a sample of red horizontal branch stars (RHB) and for this sample 



they derived the steepest metallicity gradient for the thick disc currently in the literature. 
RHB stars are very old, on average, so it is feasible that they can have a different metallicity 
variation than younger stars. However, the difference between their metallicity gradient and 
others in the literature for the thick disc is large, which motivates confirmation by other 
works with different data. 



In our previous study (jCoskunoglu et al.ll2011al . paper II), we investigated the metallicity 



gradient of a dwarf sample and we confirmed the radial metallicity gradient of - 0.04 dex 



kpc 1 based on calibrated metallicities from RAVE DR3 ( jCoskunoglu et al 



2011aL and the 



references therein, Paper I). Additionally, we showed that the radial metallicity gradient is 
steeper for our sample which is statistically selected to favour younger stars, i.e. F-type stars 
with orbital eccentricities ey ^ 0.04, i.e. d[M/ H]/dR m = —0.051 ±0.005 dex kpc -1 . Vertical 
metallicity gradients could not be derived in this study due to short z-distances from the 
Galactic plane. Therefore, in this present study, we analyses stellar abundance gradients of 
red clump stars from RAVE data. 

Red Clump (RC) stars are core helium-burning stars, in an identical evolutionary phase 
to those which make up the horizontal branch in globular clusters. However, in intermediate- 
and higher-metallicity systems only the red end of the distribution is seen, forming a clump 
of stars in the colour-magnitude diagram. In recent years much work has been devoted to 
studying the suitability of RC stars for application as a distance indicator. Their absolute 
magnitude in the optical ranges fro m My = +0.70 mag for tho se of spectral type G8 III 
to My = +1.0 mag for type K2 III ( Keenan fc Barnbauml 1999 ). The absolute magnitude 
of these sta rs in the K band is = —1.61 ± 0.03 mag with negligible dependence on 
metallicity ( Alves I 2000 ). but with a real dispersion (see below), and in the / band Mj = 
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—0.23 ± 0.03 mag, again without dependence on metallicity (IPaczynski fc Stanekl Il998l ). 
RAVE DR3 red clump giants stars occupy a relatively larger ^-distance interval than do 
RAVE DR3 dwarfs. Hence, we should have sufficient data to be able to derive abundance 
gradients for both directions, radial and vertical. As the RHB stars are on t he extended 



Chen et al 



branch of RC stars, we anticipate the results to be similar to RHB stars claimed in 
(120111 ). and so are able to test their result. 

The structure of the paper is: Data selection is described in Section 2; calculated space 
velocities and orbits of sample stars is described in Section 3. Population analysis and results 
are given in Section 4 and Section 5, respectively. Finally, a discussion and conclusion are 
presented in Section 6. 



2 DATA 



Siebert et al 



20111 ) 



The data were selected from the third data release of RAVE (DR3; 
RAVE DR3 reports 83072 radial velocity measurements for stars 9 ^ I ^ 12 mag. This 
release also pro vides stellar atmosp heric parameters for 41 672 spectra representing 39833 



individual stars ( iSiebert et al. 



20111 ). The accuracy of the radial velocities is high, marginally 



improved with DR3: the distribution of internal errors in radial velocities has a mode of 
0.8 km s _1 and a median of 1.2 km s _1 , while 95 per cent of the sample has an internal 
error smaller than 5 km s -1 . The uncertainties for the stellar atmospheric parameters are 
typically: 250 K for effective temperature T e ff, 0.43 dex for surface gravity logg and 0.2 dex 
for [M/H]. While RAVE supports a variety of chemical abundance scales, we use here just 
the public DR3 values. Since anticipated gradients are small, this provides a well-defined 
set of parameters for analysis. The proper motions of the stars were taken from RAVE 
DR3 values, which were compiled from PPMX, Tycho-2, SSS and UCAC2 catalogs. The 
distribution of RAVE DR3 stars in the Equatorial and Galactic coordinate planes are shown 
in Fig. 1. 

We applied the following cons traints to obtain a homogeneous RC sample with best 
quality: f 
(2MASS; 



2 ^ log q {cm s 2 ) ^ 3 (IPuzeras et al.ll2010l ). ii) the Two Micron All Sky Survey 



Skrutskie et al 



2006)_photometric data are of quality labelled as "AAA" , and iii) 



(J — H)q > 0.4 ( iBilir et al.ll201ll ). The numerical values for \ 2 an d the median of S/N value 
of sample spectra (totally 7985 stars), thus obtained are 671 and 41, respectively. Proper 
motions for 139 out of the 7985 stars are not available in the RAVE DR3, hence these values 
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Figure 1. Distribution of RAVE DR3 stars in the Equatorial (top) and Galactic (bottom) coordinate planes. 



are provided from the PPMXL Catalog of lRoeser. Demleitner fc Schilbachl ( 120101 ) . Distances 
were obtained by combining the apparent K s magnitude of the star in query and th e absolute 
magnitude Mk s = — 1.54±0.04 mag, adopted for all RC stars (lGroenewegenll2008l ) . Whereas 
the E(B — V) reddening were obtained iterativel y by using published met hodology (for more 



Coskunoglu et al. 



2011b , and the references 



detailed information regarding the iterations see 
therein;). Then, the de-reddening of magnitud es and colours in 2MASS w ere carried out by 
the following equations with the co-efficient of iFiorucci fc Munaril (120031 ): 



J Q = J- 0.887 x E(B - V) 
( J - H) = ( J - H) - 0.322 x E(B - V) (1) 
(H - K s ) = (H - K s ) - 0.183 x E(B - V) 

Note the real dispersion in absolute magnitude among RC stars (Fig. 2). This will have 
the effect of blurring the derived distances, and so smoothing any derived gradient. Given 
that we search for a linear gradient, such distance uncertainties will tend to somewhat reduce 
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Figure 2. The colour-abso lute magnitude diagram for the revised Hipparcos catalogue (van Lccuwcn 2007|) cross-matched with 
2MASS l|Cutri et alj|2003h . The very well illustrated red clump justifies the colour selection used in our study. 

any measured gradient. Given our results below, we consider any such effect to be second 
order. 

The distance range of the sample and the median of the distances are 0.2 ^ d ^ 3.4 
kpc and 1.34 kpc, respectively (Fig. 3), which are sufficient to investigate both radial and 
vertical metallicity gradients. The distribution of colour excess E(B — V) is given in three 
categories, i.e. 0° < \b\ ^ 30°, 30° < \b\ ^ 60°, and 60° < \b\ < 90°, whose mean values are 
0.14, 0.06, and 0.02 mag, respectively (Fig. 4). The reddening is rather small at intermediate 
and high Galactic latitudes, as expected. The projection of the sample stars onto the (X, 
Y) and (X, Z) planes (Fig. 5) show that their distribution is biassed (by design: RAVE does 
not observe towards the Galactic bulge and so), the median of the heliocentric coordinates 
are X = 0.52, Y = -0.59, Z = -0.41 kpc. 




3 SPACE VELOCITIES AND ORBITS 



We combined the distances estimated in Section 2 with RAVE kinematics and the available 
prope r motions, applying the algorithms and the transformation matrices of L 



Johnson &: Soderblom 



(119871 ) to obtain their Galactic space velocity components (U, V, W). In the calculations 
epoch J2000 was adopted as described in the I nternation al Celestial Reference System 
(ICRS) of the Hipparcos and Tycho-2 Catalogues ( )ESAlll997l ). The transformation matrices 
use the notation of a right handed system. Hence, U, V, and W are the components of a 
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Figure 3. Cumulative (panel a) and frequency (panel b) distributions of distances of our sample of RAVE red clump giants. 
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Figure 4. Distribution of colour excess E(B — V) for our sample of RAVE red clump giants. 

velocity vector of a star with respect to the Sun, where U is positive towards the Galactic 
centre (I = 0°, b = 0°), V is positive in the direction of Galactic rotation (/ = 90°, b = 0°) 
and W is positive towards the North Galactic Pole (b = 90°). 

We adopted the value of the rotation speed of the Sun as 222.5 km s -1 . Correction for 




differential Galactic rotation is necessary for accurate determination of U, V, and W velocity 
components. The effect is proportional to the projection of the distance to the stars onto 



the Gala ctic plane, i.e. the 



rotation (IMihalas fc Binney 



V ve locity component is not affecte d by Galactic differentia l 



19811 ). We applied the procedure of IMihalas fc Binneyl (119811 ) 



to the distribution of the sample stars and estimated the first order Galactic differential 
rotation corrections for U and V velocity components of the sample stars. The range of 
these corrections is —92 < dU < 58 and 6 < dV < 7 km s -1 for U and V, respectively. 
As expected, U is affected more than the V component. Also, the high values for the U 
component show that corrections for differential Galactic rotation can not be ignored. One 



notices that Ga 
ones for dwarfs 



actic differential rotation corrections are rather larger than the corresponding 



Coskunoglu et al.l ( 1201 lbf ). The U, V, and W velocities were reduced to LSR 
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by adopting the solar LSR velocities in ( ICoskunoglu et al.ll2011bl ). i.e. (U Q , V@, W Q ) = (8.83, 
14.19, 6.57) km s _1 . 

The uncertainty of the space velocity components U err , V err and W err were computed by 
propagating the uncert ainties of the proper motions , distances and radial velocities, again 
using an algorithm by iJohnson fc Soderbloml (119871 ). Then, the error for the total space 
motion of a star follows from the equation: 



"7 



2 

err 



u 2 +v' 2 +w 



2 

err 



r2 
err' 



(2) 



The mean S err and standard deviation (s) for space velocity errors are S err = 39 km s _1 and 
s = 36 km s -1 , respectively. We now remove the most discrepant data from the analysis, 
knowing that outliers in a survey such as this will preferentially include stars which are 
systematically mis-analysed binaries, etc. Astrophysical parameters for such stars are also 
likely to be relatively unreliable. Thus, we omit stars with errors that deviate by more than 
the sum of the standard error and the standard deviation, i.e. S err > 75 km s _1 . This 
removes 1204 stars, ~15.1 per cent of the sample. Thus, our sample was reduced to 6781 
stars, those with more robust space velocity components. After applying this constraint, the 
mean values and the standard deviations for the velocity components were reduced to (U err , 
Verr, W en .)=(15.03 ± 10.61, 15.12 ± 11.20, 15.68 ± 12.06) km s _1 . The distribution of the 
errors for the space velocity components is given in Fig. 6. In this study, we used only the 
sub-sample of stars (6781 stars) with standard error S err ^ 75 km s -1 . The U, V, W velocity 
diagrams for these stars are shown in Fig. 7. The centre of the distributions deviate from 
the zero points of the U, V, and W velocity components, further indicating the need for a 
Local Standard of Rest reduction. 

To complement the chemical abundance data, accurate kinematic data have been ob- 
tained and used to calculate individual Galactic orbital parameters for all stars. The shape 
of the stellar orbit is a proxy, through the age- velocity dispersion relation, for age, with more 
circular orbits hosting statistically younger stars. 



In order to calculate t 
described in the literature ( 


rose parameters we 


usee 


standard gravitational potentials well- 


Mivamoto & Naerai 


1975 


Hernauist 


1990 


Johnston. Soereel h Hernauist 


1995: 


Dinescu. Girard & van Altena 


1999 


) to estimate orbital elements of each of the sample 



stars. The orbital elements for a star used in our work are the mean of the corresponding 
orbital elements calculated over 15 orbital periods of that specific star. The orbital integra- 
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Figure 6. Error histograms for space velocity (panel a) and its components (panels b-d) for our sample of RAVE red clump 
stars. The arrow in panel (a) indicates the upper limit of the total error adopted in this work. The shaded part of the histogram 
indicates the error for different velocity components of stars after removing the stars with large space velocity errors. 



tion typically corresponds to 3 Gyr, and is sufficient to evaluate the orbital elements of solar 
neighbourhood stars, most of which have orbital periods below 250 Myr. 



Solar neighbourhood velocity space inclu des well-estab 



ble c l assic moving groups or stellar streams 



1999 



Nordstrom et al. 



200J). 



Famaey et al. 



Dehnen 



1998 



ished substructures that resem- 



Sku lian. Hearnshaw fc Cottrell 



( 120051 ) , iFamaey. Siebert fc Jorissenl (120081 ) and 
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Figure 7. The distribution of velocity components of our final cleaned sample of RAVE red clump stars with high-quality 
data, in two projections: (a) U — V and (b) W — V. 




Pompeia et al.l (120111 ) show that, although these streams include clusters, after which they 
are named, and evaporated remnants from these clusters, the majority of stars in these 
streams are not coeval but include stars of different ages, not necessarily born in the same 
place nor at the same time. They argue these streams are dynamic al (resonant) in origin, 



proba bly related to dynamical perturbations by transient spiral waves (IDe Simone. Wu fc Tremaine 



20041 ). which migrate stars to specific regions of the [/V-plane. Stars in a dynamical stream 
just share a common velocity vector at this particular epoch. These authors further point 
out the obvious and important point that dynamical streams are kinematically young and so 
integrating backwards in a smooth stationary axisymmetric potential the orbits of the stars 
belonging to these streams is non-physical. RAVE stars are selected to avoid the Galactic 
plane (|6| > 10°) . Dynamical perturbations by transient spiral waves are strongest closest to 
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the Galactic plane so there will be fewer dynamical stream stars in our RAVE sample. Hence, 
contamination of the dynamical stream stars is unlikely to affect our statistical results. 

To determine a possible orbit, we first perform test-particle integration in a Milky Way 
potential which consists of a logarithmic halo of the form 

$ halo (r)=t; 2 ln(l + ^), (3) 

with vq = 186 km s _1 and d — 12 kpc. The disc is represented by a Miyamoto-Nagai 
potential: 

<W^, z) = ° Md (4) 

R 2 + (a d + Jz 2 + b\ 



with Md = 10 11 M Q , a& = 6.5 kpc and = 0.26 kpc. Finally, the bulge is modelled as a 
Hernquist potential 

$bui gc (r) = -^J, (5) 

using Mb = 3.4 x 10 10 M Q and c = 0.7 kpc. The superposition of these components gives 
quite a good representation of the Milky Way. The circular speed at the solar radius is ~ 220 
km s -1 . Plsr — 2.18 x 10 8 years is the orbital period of the LSR and V c = 222.5 km s _1 
denotes the circular rotational velocity at the solar Galactocentric distance, Rq = 8 kpc. 

For our analysis of gradients, we are interested in the mean radia l Galactocentr i c dist ance 
(R m ) as a function of the stellar population and the orbital shape. [Wilson et al.l (120111 ) has 
analysed the radial orbital eccentricities of a RAVE sample of thick disc stars, to test thick 
disc formation models. Here we focus on possible local gradients, so instead consider the 
vertical orbital eccentricity, ey. R m is defined as the arithmetic mean of the final perigalactic 
(R p ) and apogalactic (R a ) distances, and z max and z m i n are the final maximum and minimum 
distances, respectively, to the Galactic plane. Whereas e v is defined as follows: 



( | Z m ax | | ^min \ ) 
Rm. 



(6) 



where R m = (R a + R p )/2 (|Paulill2005l ). Due to z-excursions R p and R a can vary, however 



this variation is not more than 5 per cent. 
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4 POPULATION ANALYSIS 

4.1 Classification using space motions 



The procedure of iBensby. Feltzing fc Lundstroml (120031 ) : iBensby et al.l (120051 ) was used to 



separate sample stars into different populations. This kinematic methodology assumes that 
Galactic space velocities for the thin disc (-D), thick disc {TD), and stellar halo (H) with 
respect to the LSR have Gaussian distributions as follows: 

(VlSR — Vasym) 2 



f(U, V, W) — k x exp 
where 
k 



u 2 

U LSR 



?rr 2 



(7) 



1 



L S R ® V L S R ®W L S R 

normalizes the expression. For consistency with other analyses <Ju L3B ■ 
adopted as the characteristic velocity dispersions: 35, 20 and 16 km s -1 for thin disc {D); 



(8) 

a V LS R and VWlsr were 



67, 38 and 35 


cm s 


(Bensbv et al. 


2003) 



20031 ). Vasym is the asymmetric drift: -15, -46 and -220 km s 1 for t 



l in disc , 
fhoilri ) 



thick disc and halo, respectively. LSR velocities were taken from ICoskunoglu et al. 
and these values are (U, V, W) LSR = (8.83 ± 0.24, 14.19 ± 0.34, 6.57 ± 0.21) km s _1 . 

The probability of a star of being "a member" of a given population with respect to a 
second population is defined as the ratio of the f(U, V, W) distribution functions times the 
ratio of the local space densities for two populations. Thus, 

TD/D = ^IR x flR TD/H =^IR x llR y (9) 

Xd Jd Xh /fl- 

are the probabilities for a star being classified as a thick disc star relative to it being a 

thin disc star, and relative to it being a halo star, respectively. X D , X TD and X H are the 

l ocal space densities for thin disc, thick disc and halo, i.e. 0.94, 0.06, and 0.001 5, respectively 



Robin et al. 



1996 



Buser. Rong fc Karaalilll999l ) . We followed the argument of 



Bensbv et al. 



(120051 ) and separated the sample stars into four categories: TD/D ^0.1 (high probability 
thin disc stars), 0.1 < TD/D ^ 1 (low probability thin disc stars), 1 < TD/D ^ 10 
(low probability thick disc stars) and TD/D > 10 (high probability thick disc stars). Fig. 
8 sho ws the U — V and W — V diagrams as a function of population types defined by 
using IBensby et al.l (120031 ) 's criteria. It is evident from Fig. 8 that the kinematic population 
assignments are strongly affected by space-motion uncertainties. 3385 and 1151 stars of the 
sample were classified as high and low probability thin disc stars, respectively, whereas 646 
and 1599 stars are low and high probability thick disc stars (Table 1). The relative number 
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Figure 8. U — V and W — V diagrams for our sample of red clump stars, applying Bcnsby et al. (2003)'s population classification 
criterion. It is apparent that space-motion uncertainities remain significant, even for this sample. 



Table 1 . The space velocity component ranges for the range of population types into which our RAVE red clump stars have 
been classified. 



Parameters 


U 


V 


W 


N 




(km s _1 ) 


(km s" 1 ) 


(km s" 1 ) 




TD/D s: 0.1 


(-98, 97) 


(-60, 51) 


(-44, 43) 


3385 


0.1 < TD/D ^ 1 


(-133, 123) 


(-78, 72) 


(-57, 58) 


1151 


1 < TD/D ^ 10 


(-151, 154) 


(-93, 85) 


(-69, 69) 


646 


TD/D > 10 


(-337, 290) 


(-400, 291) 


(-344, 273) 


1599 



of high probability thick disc (RC) stars are much larger than the corresponding ones in 
Coskunoglu et al] ( 2011a ) (2 per cent), i.e. 24 per cent. 



4.2 Population classification using stellar vertical orbital shape 

Both radial and vertical orbital eccentricities contain valuable information: here we con- 
sider the vertical orbit shape. Vertical orbital eccentricities were calculated, as described 
above, from numerically- integrated orbit s. We term th i s the dynamical method of popula- 

(120031 ) 's approach. The distribution 



tion assignment, which complements the 



Bensby et al. 



function of ey is not consistent with a single Gaussian distribution, as is shown in Fig. 9. A 
two-Gaussian model however does provide an acceptable fit. For convenience, we separated 
our sample into three categories, i.e. stars with e v ^ 0.12 (3448 stars), 0.12 < e v ^ 0.25 
(2389 stars) and e v > 0.25 (944 stars), and fitted their metallicities to their mean radial dis- 
tances (R m ) in order to investigate the presence of a metallicity gradient for RAVE RC stars. 
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Figure 9. The distribution of vertical orbital eccentricity for our sample of RC stars. The dashed lines indicate to the stars 
with ey ^ 0.12 (the one with higher mod) and 0.12 < ey ^ 0.25 eccentricities (lower mod), respectively, whereas the solid line 
corresponds to the whole sample. 

We provide in Table 2 (electronically) for each star, stellar parameters from RAVE DR3, 



calculated kinematical and dynamical parameters and our stellar population assignment. 



Table 2. Stellar atmospheric parameters, astrometric, kinematic and dynamic data for the whole sample: (1): Our catalogue number, (2): RAVEID, (3-4): Equatorial coordinates in CTi 
degrees (J2000), (5): T e g- in K, (6) logg (cm s~ 2 ) in dcx, (7): Calibrated metallicity [M/H] (dcx), (8-9): Proper motion components in mas/yr, (10): d in kpc, (11): Heliocentric radial 
velocity in km s -1 , (12-17): Galactic space velocity components, and their respective errors in km s — 1 , (18-19): Perigalactic and apogalactic distance in kpc, (20-21): Minimum and ^ 
maximum distance from the Galactic plane in kpc, (22): TD/D ratio as mentioned in the text. TD/D = 999.999 shows the value of TD/D > 999. 
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Figure 10. The metallicity distribution for our cleaned sample of RC stars, using RAVE DR3 abundances. 
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Figure 11. The normalized metallicity distribution as a function of assigned (probabilistic) population types. 

5 RESULTS 

5.1 First hints of a metallicity gradient apparent in the RC sample of stars. 

We show in Fig. 10 the distribution of metallicities for our final sample of RAVE RC stars. 
The metallicity distribution for all populations is rather symmetric with a mode at [M / H] ~ 
—0.4 dex. In Fig. 11 we show the normalized metallicity distribution functions, with the 
sample sub-divided by probabilistic population assignment, as described above. This Fig. 11 
gives an indication of a systematic shift of the mode, shifting to low metallicities when one 
goes from the thin disc stars to the thick disc stars (Fig. 11). The ^-distance distribution of 
our sample is shown in Fig. 12. While the typical star is at distance \z\ ~ 0.5 kpc, there is a 
significant sample at larger distances. The range in z— distance is large enough to allow us 
to consider vertical metallicity gradient estimation. 

The whole sample of stars, with no consideration of population assignment, were sepa- 
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Figure 12. Distribution of the distances from the Galactic plane for our red clump sample 



Table 3. The mean metallicity and vertical eccentricities as a function of distance from the Galactic plane. 



z ranges 


N 


< z > 


< M/H > 


< e v > 


(kpc) 




(kpc) 


(dex) 




(0.0, 0.5] 


2261 


0.35 


-0.33±0.25 


0.087±0.106 


(0.5, 1.0] 


3376 


0.70 


-0.41±0.25 


0.161±0.136 


(1.0, 1.5] 
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-0.47±0.25 
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1.71 
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0.378±0.222 


(2.0, 2.5] 


54 


2.20 


-0.56±0.24 
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rated into five bins in distance, and mean z-distances from the Galactic plane, mean metal- 
licities, and mean vertical eccentricities were calculated for each bin. These are presented 
in Table 3. Fig. 13 summarises the dependence of the results on vertical distance from the 
Galactic plane. The apparent variation of the mean metallicity with ^-distance from the 
Galactic plane in Fig. 13 indicates the existence of a vertical metallicity gradient for RC 
stars. This figure does not however allow discrimination between a true gradient in a consis- 
tently defined population (e.g. "thick disc") and a changing relative contribution from two 
populations (e.g. "thin disc" and "thick disc"), which have different modal abundances, and 
for which either one, both, or neither has an intrinsic gradient. Also given in Fig. 13 is the 
variation of the vertical eccentricity with z, which again shows either a smooth transition 
from thin disc to thick disc eccentricities, or a changing population mix, or both disc. 



5.2 Metallicity gradients using the kinematical population assignment method 

We consider the metallicities as a function of the mean orbital Galactocentric radial distance 
(R m ) and maximum distance to the Galactic plane (|-2 ma a;|) for each different population 




defined in Section 4, and test for the presence of vertical and radial metallicity gradients 
for each population. We fitted the distributions to linear equations, whose gradient is any 
metallicity gradient, d[M/H]/dR m or d[M/H]/dz max . The results are shown in Fig. 14. 

The radial d[M/H]/dR m gradients are small or consistent with zero. The best deter- 
mined value (largest ratio of gradient value to error value) is for low probability thin disc 
stars, where the gradient is d[M/H]/dR m = —0.031 ± 0.003 dex kpc -1 (Table 4). The only 
metallicity gradient consistent with zero (with small formal errors) is for low probability 
thick disc stars. 

The vertical metallicity gradients are (absolutely) much larger than the radial ones, 
and statistically are detected. The largest gradient in the vertical direction is steeper for 
high probability thin disc stars, relative to the other populations, viz d[M/H]/dz max = 
—0.109 ± 0.008 dex kpc -1 . Additionally, the metallicity gradient for high probability thick 
disc stars is not zero, i.e. d[M/H]/dz max = —0.034 ± 0.003 dex kpc -1 . 
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Table 4. Radial and vertical metallicity gradients for red clump RAVE stars evaluated from kinematical and dynamical data. 
The meanings of the probabilistic population assignments, TD/D and e v , are explained in the text. 



Population type 


d[M/H]/dR m 


d[M/H]/dz max 


Sample size 




(dex kpc ) 


(dex kpc ) 




TD/D sC 0.1 


-0.041±0.003 


-0.109±0.008 


3385 


0.1 < TD/D < 1 


-0.031±0.003 


-0.086±0.013 


1151 


1 < TD/D < 10 


-0.001±0.005 


-0.036±0.016 


646 


TD/D > 10 


0.017±0.008 


-0.034±0.003 


1599 


e v < 0.07 


-0.041±0.007 


-0.260±0.031 


1269 


e v < 0.12 


-0.025±0.040 


-0.167±0.011 


3448 


0.12 < e v s: 0.25 


-0.013±0.004 


-0.103±0.008 


2389 


e v > 0.25 


0.022±0.006 


-0.022±0.005 


944 



5.3 Metallicity gradients using the dynamical population assignment method 

We now consider the metallicities as a function of the mean orbital Galactocentric radial 
distance (R m ) and maximum distance to the Galactic plane ( | ^maa; | ) for the populations 
defined by their eccentricities in Section 4.2, with one slight modification, in that we add 
a population defined by e v ^ 0.07 (1269 stars) representing the blue st ars. We make this 



modific ation due to our experience in analysing the RAVE dwarf stars. In ICoskunoglu et al. 



(j2011al ). we showed that the blue stars which have the smallest orbital eccentricities (most 
circular orbits) have also steeper metallicity gradients relative to samples with larger orbital 
eccentricities. 

The results are presented in Table 4 and Fig. 15. The steepest metallicity gradient is 
for e v ^ 0.07 in the vertical direction, i.e. d[M/H}/dz max = —0.260 ± 0.031 dex kpc -1 . The 
vertical metallicity gradient systematically decreases with increasing e v values, becoming 
close to zero for the largest vertical eccentricities, e v > 0.25. As noted above, our (red 
clump giants) sample consists of thin and thick disc stars. The largest orbital eccentricities 
correspond to thick disc stars. 

That is, the vertical metallicity gradient for the thick disc is close to zero. The trend of 
the radial metallicity gradient is almost the same as the vertical metallicity gradient; less 
steep. For example, d[M/H]/dR m = -0.041 ± 0.007 dex kpc -1 for e v ^ 0.07. 



6 DISCUSSION AND CONCLUSION 

We have used the RAVE DR3 data release to identify RC stars, further excluding cool 
stars, and those with the most uncertain space motions. We used the calibrated RAVE DR3 
metallicities and the mean radial and maximum distances to investigate the presence of 
radial and vertical metallicity gradients, dividing the sample into a variety of subsamples. 
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Figure 14. Radial (R m -[M / H]) and vertical (z m ax-[M / H]) metallicity gradients for the red clump sub-samples, subdivided 
into probabilistic population types as described in the text. 



We derive significant radial and vertical metallicity gradients for high probability thin 
disc stars and for the subsample with e v ^ 0.07. We derive significant and marginally 
shallower gradients for the other subsamples. We do not detect any significant gradients 
for thick disc stars. Vertical metallicity gradients are much steeper than the radial ones, 
for the same subsample, as (perhaps) expected. We derive the metallicity gradients for the 
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Figure 15. Radial (R m -[M//f]) and vertical (z ma x-[M / H]) metallicity gradients for the red clump sub-samples, subdivided 
by vertical orbital eccentricity. 



subsample e v ^ 0.07 which are the youngest sample stars, to be d[M/H]/dR m = —0.041 ± 
0.007 dex kpc" 1 and d[M/H]/dz max = -0.260 ± 0.031 dex kpc" 1 . 

The radial metallicity gradients for all subsamples are rather close to the corresp onding 



ones obtained for F-type dwarfs. Hence, the discussion in 



Coskunoglu et al. 



(2011a 



) is also 



valid here. We cannot determine any vertical metallicity gradient for RAVE dwarfs due to 



their small distances fro m the Galactic p 



those obtained for giants. 



Chen et al. 
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ane. However, we can compare our results with 
(1201 ll ) investigated the metalli city gradient of the thick 
disc by using Red Horizontal Branch (RHB) stars from SDSS DR8 (jAihara et al.ll201l[ ) and 
they found two different vertical metallicity gradients estimated in two ways. One is a fit 
to the Gaussian peaks of the metallicity histograms of the thick disc by subtracting minor 
contributions from the thin disc and the inner halo based on the Besangon Galaxy model 
(IRobin et al.lll996f ). The resulting gradient is d[M/H)/dz = -0.12 ± 0.01 dex kpc -1 for 
0.5 < \z\ < 3 kpc. The other method is fitting the data linearly for the stars 1 < \z\ < 3 
kpc where the thick disc is the dominant population. Five subgroups were then selected in 
different directions in the X — \z\ plane to investigate the difference in the vertical metallicity 
gradient between the Galactocentric and anti-Galactocentric directions. They found that a 
vertical gradient of d[M/H]/dz = —0.22 ± 0.07 dex kpc -1 was detected for five directions 
except for one involving the population of stars from the b ulge. 



Chen et al 



(120111 ) are in agree- 



Neither of the vertical metallicity gradients claimed by 
ment with our results, nor are they consistent with studies appearing in previous works in 
the literature. 



of 



Another rec e nt inv estigation of the vertical metallicity gradient for the thick disc is that 



Ruchti et al. 



(120111 ). In that paper the authors gave a sample of 214 Red Giant Branch 
stars, 31 red clump/red horizontal branch stars, and 74 main- sequence/ sub-giant branch 
metal-poor stars. They found that the thick disc [a/Fe] ratios are enhanced, and have little 
variation (< 0.1 dex). Their sample further allowed, for the first time, investigation of the 
gradients in the metal-poor thick disc. For stars with [Fe/H] < —1.2 dex, the thick disc 
shows very small gradients, < 0.03 ± 0.02 dex kpc -1 , in a enhancement, while they found a 
d[Fe/H]/dR = +0.01±0.04 dex kpc -1 radial gradient and a d[Fe/H]/dz = -0.09±0.05 dex 
kpc -1 vertical gradient in iron abundance. We consider only the gradient in iron abundance, 



not the gradient in a enhancement. We may transform pub 



metallicity values by means of the equation ( iZwitter et al. 

[M/H] = [Fe/H] + 0.11 [1 ± (1 - e -3.6|[Fe/H]+0.55|^ 



ished iron abundances to RAVE 



2008( ) 



(10) 



This reveals that the 



Ruchti et al. 



(120111 ) vertical gradient is consistent with the vertical 
metallicity gradient determined here, within the errors, for high probability thick disc stars, 
d[Fe/H]/dz max = —0.034 ± 0.003 dex kpc -1 . However there is a difference between the 
corresponding radial metallicity gradients in the two studies. An explanation for this dis- 
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agreement may be the differing metallicity range of the sample stars used in the two works. 
In the present study we consider stars with [M/H] > —1.1 dex, with a minority at the 
metal-poor tail in our work, while the corresponding selection is [Fe/H] < —1.2 dex in 



Ruchti et all (120111 ). 

The radial metallicity gradients we have estimated from red clumps giants stars are 



consistent with those derived in paper II for dwarfs (ICoskunoglu et al.l l2011af ). The ro- 
bust metallicity gradients we determine are d[M / H]/dR m = —0.041 ± 0.003 dex kpc -1 for 
the high probability thin disc stars, the population type labelled with TD/D ^ 0.1, and 
d[M/H]/dR m = —0.041 ± 0.07 dex kpc -1 for the sample with eccentricity e v ^ 0.07. Sam- 
ples biased to low probability thin disc and thick disc stars show systematically shallower 
gradients. Complementary to the dwarf sample, the distance range of the red clump giant 
stars, with median distance of 1.4 kpc, provides information on vertical metallicity gradi- 
ents. The vertical metallicity gradients for the high probability thin disc stars and for the 
sample with e v ^ 0.07 are d[M/H}/dz max = -0.109 ±0.008 dex kpc" 1 and d[M/H}/dz max = 
—0.260 ±0.031 dex kpc -1 , respectively. For high probability thick disc stars we could detect 
a vertical metallicity gradient; a shallow one however, i.e. d[M/H]/dz max = —0.034 ± 0.003 
dex kpc -1 . 

From our analysis, we may conclude that, despite their greater distances from the Galac- 
tic plane, the RAVE DR3 red clump giant stars confirm the radial metallicity gradients 
found for RAVE DR3 dwarf stars. Because of their greater distances from the Galactic 
plane, the RAVE DR3 RC stars also permit vertical metallicity gradients to be measured. 
These findings can be used to constrain formation scenarios of the thick and thin discs. 
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